Screening of the Photon Field in Surface Photoemission from 

Simple Metals 

S. Natschlager^, J. Alberich^, N. Barberan^, J. Bausells^, R. Ruppin^ and C Simo^ 

^ Institut fur Theoretische Physik, Johannes Kepler Universitdt, A-404O Linz, Austria 
f^ , "^Departament d'Estructura i Constituents de la Materia, Facultat de Fisica, Universitat de 

^^ ! Barcelona, E- 08028 Barcelona 

^Centro Nacional de Microelectronica-CSIC, E-08193 Bellaterra (Barcelona), Spain 
(V-i : '^SOREQ NRC, Yavne 81800, Israel 



(N 
> 



^ Departament de Matemdtica Aplicada i Analisi, Universitat de Barcelona, E-08007 

Barcelona, Spain 



o 

g;j ■ Abstract 

-)— > 

^ [ A previous two-step hydrodynamic model calculation of angle-resolved 

'^ I photoemission has been improved by modeling the surface with a soft ana- 

Q I lytic function. The same ground-state potential is used to screen the photon 

- ^_, , field and to calculate the electronic wave functions. The implemented code 

H , can be used for arbitrary ground state densities along the perpendicular 

direction. Mathematical and physical aspects of the model are discussed 
extensively. The behavior of the main peak of the spectrum in response to 
the variation of some parameters leads us to the conclusion that the peak 
cannot be interpreted as the direct excitation of multipole modes. This 
conclusion was suggested in the previous two-step model calculation. 
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1. Introduction 

The subject of the interaction of external probes with metal surfaces, and the pos- 
sibility of excitation of multipole modes, have attracted much attention. The multipole 
modes, which were first predicted theoretically by BennettGJ in 1970, from the hydro- 
dynamic model, have only recently been detected experimentally by inelastic electron 
scattering on Na, K, Cs and Ag surfacesfla Theoretical quantum-mechanical interpreta- 
tions of these experiments, based on the random phase approximation and on the local 
density approximation, have been presented.l3i3'El Although no multipole modes could be 
observed in electron scattering from Al surfaces, it was debated whether a peak which ap- 
pears in the photoemission spectrum of AlQ'Ll could be attributed to multipole modes.&Ero 
To provide more insight into this question, theoretical calculations of the angle depen- 
dence of the photoemission spectrum are required. Here we present such a calculation, 
using a refined hydrodynamic model, into which a realistic surface density profile is 
incorporated. 

Photoemission within a solid needs some linear momentum supplier besides the two 
components: the incoming photon and the emitted electron. An isolated electron cannot 
absorb a photon in a process that preserves total energy and momentum. Due to the fact 
that the electrons have a finite mean free path, the photon is absorbed near the surface. 
Within this region, the most important momentum suppliers in a simple metal are the 
step potential and the screened photon field. 

These two contributions can be determined explicitly in a Hamiltonian formulation. 
The perturbation produced by the incident photon is given by 

where only terms up to first order in the vector potential A are considered and where 
the gauge $ = ($ is the scalar potential) has been used, p is the linear momentum 
operator and c, e, and m are the velocity of light in vacuum and the electron charge and 
the mass, respectively. We will employ atomic units throughout the paper. 



If I z ) and I / ) are two eigenstates of the unperturbed Hamiltonian Hq and A is 
the external plus the induced field, the amplitude of the excitation probability is given 
(apart from constant factors) by: 

(/ I p ■ A + A ■ p I i ) = -i (f I (V ■ A) + 2A ■ V I i ). (2) 

While the first term on the right hand side is directly related to the spatial variation 
of the vector potential due to screening effects within the electron gas, the second term 
is related to the step potential at the vacuum-solid interface. This follows from the fact 
that for the unperturbed Hamiltonian Hq = p^ /2m + V we have 

(/|V|i) = -J--(f|VV|i). (3) 

Therefore, if A is taken out before the integral in the second term of the right hand side 
of Eq. (^ (as in some approximations with optical photons), then the entire contribution 
to the momentum would come from the variation of the ground state potential at the 
surface. However, in a full calculation, as A must be placed inside the bracket, the two 
effects (screening and inhomogeneities of V) are mixed. 

These two contributions to the excitation probability are connected with different 
kinds of processes during the electron emission: collective and single-particle excitations 
respectively. 

The contributions connected with the potential gradient term, VV, consist of indi- 
vidual electron-hole pairs created at the surface and of the subsequent emission of the 
electron. On the other hand, in the processes connected with the divergence of the vec- 
tor potential the electron gas screens the incoming photon, producing surface collective 
excitations that supply the extra momentum. 

Moreover, different types of collective excitations play different roles in photoemission. 
When the photon frequency is equal to tUp, the bulk plasma frequency, the amplitude 
of the vector potential is drastically reduced inside the metal and no photoelectrons are 
produced, giving a minimum in the photoemission spectrum. 

Ordinary surface plasmons of frequency uj^/ \/2 cannot be produced by an incident 
photon in a perfect surface, as the lightline and surface plasmon dispersion relations have 



no points in common on the uj vs. q plane due to retardation effects (g is the momentum 
parallel to the surface). 

Multipole modes associated with oscillating charge fluctuation in the direction per- 
pendicular to the surface are the ideal candidates to be the main components in the 
non-zero contribution related to the divergence term. 

It has been proved that, in contrast to the ordinary surface modes, multipole modes 
are very sensitive to the electronic structure at the surfacetl. 

The ordinary surface modes reflect bulk properties and at g = their energy has only 
Coulomb contribution. However, even at g = multipole modes have kinetic contribution 
in addition to the Coulomb part, giving information about local properties.ty Several ap- 
proximations of the surface region have been published in an attempt to model the most 
important part of the response of the system. Some use the hydrodynamic modelo'cJil^ 
and others use microscopic calculationgi^ta in which the photon field is screened with 
the same potential step as the one used to obtain the electronic wave functions. 

Here we improve a previous two-step profile calculation^ using a soft electronic profile 
to avoid the discontinuities at the vacuum-layer and layer-bulk boundaries that could be 
related with artificial enhancement of the photoemission spectrum. 

A self-consistent screening of the field is performed using an analytical expression for 
the ground state density. 

This paper is organized as follows. The photoemission formulation is presented in 
Sec. 2. Section 3 contains the calculation of the screened field and Sec. 4 contains the 
results and conclusions. 

2. Photoemission Formulation. 

We consider a semi-infinite simple metal with its surface perpendicular to the 2;-axis 
and confined to the 2; < half-space. The jellium model is assumed for the positive ionic 
density defined as 

N{z) = ho Q{-z) , (4) 



where uq = S/Anr^ is the constant bulk density, 9 is the step function and r^ is the 
electron radius. 

Taking advantage of the translational symmetry on the xy-plane, the photoemission 
flux of electrons that reach the detector per unit of the solid angle and per incident 
photon can be written as 



dl pf cos{9) 
dfl uo 



i{q + 2p\\i) j dz(t)ff(t)i + / dz(t)fi(j)i + 2 / dzcpjgcp] 
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(5) 



where we have followed the scheme developed by Mahan in Ref. |T^ (with the addition 
of the term V ■ A 7^ 0, not considered in Mahan's formulation). In Eq. (^) 0j is the 
2;-dependent part of the initial electronic state at the Fermi energy, 

$,(r) = e^P|i'-''0i(z), r=(p,z), (6) 

where pyf = pyi + q is the balance between the final and initial electronic momenta 
parallel to the surface and the parallel momentum q supplied by the incident photon. 

$/(r)=e-'Pir^0f(z) (7) 

comes from the Green function of a free electron that propagates from the metal to the 
vacuum through the surface potential at a constant energy Ef, given by 

Ef = E, + uj , (8) 

where Ei is the Fermi energy. 

The functions $j and $/ are solutions to the Schrodinger equation for the unper- 
turbed Hamiltonian Hq. 

The functions / and g are the x— and 2;— components of the screened electric field 
given by 

E = ^A, (9) 

and j = g' is the 2;-derivative of the z-component. The field components inside the 
integrals (/, g and j) are normalized to unit amplitude for the x-component of the 



incident photon in vacuum. We consider ]?-polarized light (s-polarized hght would not 
produce a longitudinal component inside the metal) with an angle of incidence 9 with 
respect to the normal to the surface. 

Two modifications to Eq.(2.9) in Ref. O have been applied to obtain our Eq.(H). 



Both of them mimic the experimental conditions used to compare our results. The first 
is that we do not perform the integration over initial states as we keep the initial state 
fixed and change 0/ as we scan the frequency. The second is the normalization of the 
output to the incident flux (photons/sec), given by 

El cos(^) 



UJ 



(10) 



where Eq is the amplitude of the incident photon far from the surface in vacuum. 

We stress that the electron crosses the surface as a free particle and the full many- 
body effects are considered in the dressed field. 

By far the most cumbersome part of the calculation is the screening of the photon 
field. We will focus on this part in the next section. 

3. Screening of the Photon Field. 

Our initial system of equations, within the hydrodynamic model, is given by 

(11) 
(12) 

(13) 

for the three unknowns E, j and n; n and j are the electronic density and current 
respectively (j = nv), A^ is the ionic density given by Eq. (^) , /5 = |f|., vf{z) being the 
local Fermi velocity and 7 = ujp{z)/d is the damping parameter. The value c? = 50 has 
been used in the present work. 

In order to linearize the equations we consider 



■E)-V'E = -^ 




V • E = 47r(n - A^) , 


# = -E- 
dt 


/3Vn - 7J , 



n = riQ + 111 exp{—iujt) , (14) 

E = Eo + El exp(-icjt) , (15) 

V = viexp(— icjt) , (16) 

l3 = l3o + l3iexp{-tujt) , (17) 

7 = 7o + 7i exp(-iu;t) , (18) 

which gives by substitution in the previous system, 

V(V • Eo) - V'Eo = , (19) 

V ■ Eo = 47r(no - A^) , (20) 

noEo = /^oVno (21) 

to zeroth order in the perturbation, and 



V ■ El = 47rni , (23) 



V(V ■ El) - V'Ei = ^Ei + ^^novi , (22) 



- 2u;noVi = noEi + hiEq - /3oVni - /9i Vno - 7onoVi (24) 

to first order in the perturbation. 

The standard procedure would be to solve the zeroth order system obtaining the 
ground state density Uq from a given ionic distribution. Instead, we proceeded in a dif- 
ferent way: we chose a well justified analytic electronic ground state profile and calculated 
Eo from Eq. (^. 

In the case of aluminum, the analytic form 

with S = 0.66 a.u. agrees well with the profile obtained from the Kohn Sham calculation 



presented in Ref. see Fig. 1. This electronic profile together with 

V{z) = -jEodz, (26) 

where 



Eo = |(37r^)^/=^ no no^/^^e, , (27) 
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obtained from Eq. (^), completely characterize the ground state of the system. In Eq. 

(p^) e^ is a unit vector along the z-axis and Uq = dno/dz. 

After some manipulations in the first-order system, we obtain the equations 

«i/ +a2f" +a^9' =0 , (28) 

bif +hf' +hg +b,g' +hg" =0 ; (29) 

the 2;-dependent coefficients a^ and 6j are given in the Appendix. 

This coupled system of second order differential equations can be transformed into 
an equivalent fourth order equation for one variable or, alternatively, into a first-order 
system of equations for the four functions yi{z) ,...,?/4(z), defined as / = yi, f = y2, 
g = ys and g' = y4, which satisfy 

y'liz) = y2iz)quad, (30) 

y'^iz) = Aiyi{z) + A^yi{z)quad, (31) 

1/3(2;) = y4iz)quad, (32) 

y'^{z) = Biyi{z) + B2y2{z) + B^ys^z) + Biyi{z)quad, (33) 

which must be solved numerically. The z— dependent coefficients Ai and Bi are defined 
as 

Ai = -ai/a2 Bi = -h/h (34) 

and are given in the Appendix. The fields were integrated from inside using the "odeint" 
and "stiff" subroutines given in Ref. ^. These codes change the size of the integration 
step according to the nature of the profile and appeared to be the appropriate method 
for our stiff profile. 

The main difficulty lies in the consideration of the boundary conditions. Deep inside 
the metal at zi < 0, where the electronic density becomes a constant function, the four 
independent analytical solutions are known: two longitudinal and two transverse plane 
wave functions moving in opposite directions. 

8 



Let us define the coefficients Ci within a general solution inside this region as 



f{z) = Cie-'"'' + C2e-""'' + Cse'"'' + C^e 



ik'iz 



^ik\z 



^ik2Z 



(35) 



where 



k2 



^(^ 



Co", 



nl/2 



pO 



uj{uj + z7o)^ 



^(^ + «7o)-WpO 2 

A ' 



1/2 



(36) 
(37) 



are z— components of the transverse and longitudinal wave vectors at Zi. The roots with 
\vii{ki) > have been used. 

The absence of incoming waves from — oo gives the ffist two of the boundary condi- 
tions, i.e., C3 = C4 = 0. 

Far from the surface in vacuum at ^2 > 0, where the density can be considered as a 
negligible constant function, the general analytical solution is given by 



ikz 



f{z) = Die'^' + 026-'"' + Dae* ^" + 0^6 



3*^3-2 



-ik^z 



(38) 



where k and k^ are the z— components of the transverse and longitudinal wave vectors 
respectively. They are obtained from equations similar to Eqs. (PB|), and (^) using the 
values of cjp, 7 and f3 at Z2. We used the condition n{z2)/no = 10^^ and zi = —Z2 to 
fix the Zi values. The remaining two boundary conditions must be the cancelation of 
both longitudinal solutions at Z2, or D^ = D^ = 0. In this way, the allowed radiation in 
vacuum is the incoming and reflected photon. 

The application of boundary conditions at different points in space does not always 
have a solution, as the system must be compatible with them. In a numerical prob- 
lem, this compatibility is a problem of relative accuracy. Taking advantage of the linear 
character of the system of equations, we inspected the 4x4 matrix M that relates the 
analytical solutions at zi and Z2 and is given by 



D = MC 



(39) 



The result of the inspection was that the matrix has a 2 x 2 box of nearly zero elements, 
those denoted by mai, 77132, m4i and 17142- This box produces negligible D3 and D4 co- 
efficients for arbitrary values of Ci and C2. However, we verified that the results were 
much more stable if we nonetheless imposed the two last boundary conditions using the 
"shooting method" .Ii3 

Figures 2 and 3 show the 2;— component of the vector potential as a function of 
z, normalized to the constant transverse component at zi. In addition to the general 
discussion that we will include in the next section, some comments about the behavior 
of the fields due to screening effects in the region where the electronic density is negligibly 
small must be included here. 

Rapid oscillations of the z— component of the vector potential appear in the region 
where the electronic density has nearly disappeared (see Fig.l for comparison). The 
main reason for this unphysical behavior is attributable to the hydrodynamic model we 
are using, which fails to reproduce the correct response when the electronic density is 
negligible. For the longitudinal component (which represents the collective excitation of 
the system), the dispersion relation in a homogeneous system is given by 

uj{uj + t-f)=ul + (3{q' + kl). (40) 

As the density goes to zero locally, /3, 7 and Up also tend to zero. However, u and q 
are finite quantities fixed by the incoming photon with which the system resonates. The 
only way to satisfy this equation is by a strong increase in k^; indeed, it must grow to 
infinity to compensate for the vanishing of the (3 coefficient. However, the consequences 
of these oscillations are not important, since the amplitude of the field tends to zero. 
This assertion is proved in the comparison shown in Fig. 4, as discussed in the next 
Section. 

This annoying effect could be avoided by the introduction of an extra term in the 
dispersion relation, the well known K^/A (K = (q, ka)) term,c2l related with electron-hole 
excitations (and coming from the kinetic term within the Hamiltonian) , which dominates 
the dispersion relation for high energy and large momentum transfer. If this term is 
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included, as the ground state density tends to zero, k^ converges to a finite value given 
by kl = 2uj — q^ and the strong oscillations would disappear. 

As was recently suggested in Ref. ^ this term can be obtained within the hydro- 
dynamic model for a homogeneous electron gas from the material equation if the local 
equilibrium approximation is relaxed and extra terms coming from higher momenta of 
the distribution function are included. We tried with two different extra terms within 
the material equation built up from the polarization vector (V ■ P = — n), one of them 
free of shear forces and given by V(V^n), and the other one containing shear forces 
and given by (V^)^P. These two terms include all the possible vector components that 
contain four derivatives of P. They both imply the dispersion relation given by: 

uj{u + ^7) =^l + /5(g' + kl) + i(g2 + klf (41) 

which reproduces the random phase approximation result of Ref. |2^. Both produce a 
new system of six coupled first-order equations which need six boundary conditions. The 
difficulty comes from the fact that two of the four solutions of k-^ in Eq. ( PD are unphys- 
ical, as for them Im(A;3)-Re(fc3)< 0, which corresponds to a plane wave that diverges in 
the forward direction producing divergent solutions as the differential equation is solved 
from Zi to Z2. The goal of our future work is to find the appropriate term in the material 
equation that would leave the system linear and simultaneously produce well-behaved 
solutions of the dispersion relation. 

4. Results and Conclusions. 

The physical parameters of our problem are the width of the surface profile given by 
5, the angle 9 between the incident photon and the normal to the surface, the angle e 
between the outcoming electron and the normal to the surface and the electron radius 

In Figs.2 and 3 the 2;— component of the vector potential is shown as a function of z 
for different values of uo. It has been normalized to the transverse component at Zi and 
consequently the deviation from 1 (of the real part) gives an idea of the importance of 
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the longitudinal component deep inside the metal. It can be observed that even though 
Im(A;2)^ Im(/i;i) and so strong damping is expected, its presence is significant at large 
distances from the surface and cannot be ignored for u > ujp. 

Some characteristic features that were previously reported by FeibelmanEj within a 
self consistent calculation are reproduced here by us: uo = uop i^ the critical value from 
which the wave vector of the longitudinal component has a nonvanishing real part and 
consequently it penetrates the metal, producing decaying oscillations towards — cxo. For 
this value of uj the amplitude of the nearly structureless field shows an abrupt decrease, 
producing a minimum (it should be zero if 7 = 0, i.e. if no damping is considered) in 
the photoemission spectrum. For uj < uOp there is a strong peak near z = produced by 
screening, a feature that strongly deviates from the classical optical calculation. Corre- 
lating the strength of the surface photoeffect with the intensity of the vector potential 
within the surface region, one would expect it to be large below ujp and much smaller 
above Up. 

Indistinguishable results were obtained when the selfconsistent density shown in Fig.l 
was used. 

In order to see the infiuence of the oscillations of the field at the vanishing density 
region on our results, we compare in Fig. 4 for 5 = 0.2 a.u. (which mimics an abrupt 
electronic profile) the evolution with u of the two transmission coefficients, Ci (transverse 
component) and C2 (longitudinal component) and the refiection coefficient Di with those 
obtained analytically from a calculation that assumes a single step, structureless surface. 
The analytical coefficients are given in Ref. |2^ by 

C,..4^f , (43) 



ek - h + (e - l)qyk2 i^ 
qt V c2 



D,^= '^-^'^\:-'''^'^\ l^-q^ , (44) 



where 



t = ek + k^-{e- l)q^/k2 , (45) 
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ujI 



^ = 1- . f . , , (46) 

and k, ki, k2 are defined in Eqs. (p6D, (p^) and (pHf). The curves do not show any structure, 
except at Up and at C(Jp/cos(^), where the transverse component can penetrate the sohd. 
The coincidence of the two results proves that the strong oscillations at the edge region 
shown in Figs. 2 and 3 do not affect the behavior of the field amplitudes as functions of 
cu. 

The same coincidence is obtained if r^ = 6 is used, keeping all the other parameters 
fixed. 

Figure 5 shows the change of the longitudinal transmission coefficient C2 [see Eq. ( p5D 
] as the electronic profile becomes less abrupt. As 6 increases, giving a fiatter profile, a 
wide peak at about u; = 0.75 Up evolves. Apparently, the position of the peak is nearly 
insensitive to the changes in the electronic profile. 

Finally, we centered on the study of the photoemission intensity per unit of the solid 
angle and incident photon for j)-polarized light as a function of the photon frequency. 
The initial state $i is at Ei = —0.1615 a.u. . This energy was chosen to reproduce the 
threshold value of u for aluminum. We take the zero of energy to be at vacuum level. 
No final state effects are taken into account. This last approximation is supported by 
experiment as no qualitative difference is observed between photoemission spectra made 
on different surfaces of aluminum crystals, suggesting that in a nearly free electron metal 
the screening effects are much greater than the structure produced by band effects.^ 

The finite mean free path of the emitted electron inside the metal has been modeled 
by a decreasing exponential factor multiplying the initial electronic wave function which 
ensures that only electrons coming from the last 50 A will be considered. 

Figure 6 shows the main result that must be compared with the experimental 
outputs^ (for 6 = Tc/A, e = and 6 = 0.66 a.u.). The main peak is centered at 
u! = 0.77 Up, in excellent agreement with the experimental peak position. The second 
peak, at uj > ujp, has a shape comparable to the experimental one, though it has greater 
intensity. This too high intensity is caused by the absence of electron-hole excitation 
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processes from the hydrodynamic model. Thus, our model lacks the Landau damping 
mechanism, through which the collective excitations can decay into electron-hole pairs. 
This type of damping would have a significant effect only for uj > Up, where the bulk 
plasmons created by the incident photon can overlap the electron-hole band and dis- 
appear, thus causing the photoemission efficiency in this region to decrease. Unlike the 
surface modes at u < Up, these bulk plasmons can have a high enough perpendicular 
momentum, to enable this process. 

More recent experimental data reported in Ref. |^ show the main peak at cj = 0.85 



Up. 



To gain insight into the physical significance of the main peak, we fixed 6 = 1.2 a.u., 
e = and Vg = 2.07 and looked for the variation of the spectra for different values of the 
direction of the incident photon (given by ^). A large value of 6 was chosen to increase 
the probability of mulitpole modes. As shown in Fig. 7, the peak position turns out to 
be sensitive to the value of 6. The peak disperses to lower energies as q\\ increases. 

A peak produced by a multipole mode should hardly disperse (for these small values of 
q) and, in any case, according to previous calculations,^ it would disperse in the opposite 
direction. Our results indicate that the main peak need not be a direct manifestation of 
multipole modes, as has been asserted in the literatureja using the coincidence in energy 
as the main argument. 

Another known difficulty concerning the identification of the main photoemission 
peak with the multipole modes, is that the RPA calculation of the photoemission 
spectrun£3 predicts that the frequency of this peak should increase with r.,. However, the 
results of measurements of the multipole modes frequency by electron loss spectroscopy 
do not follow this trend. 

In order to understand the origin of the 6 dependence of dl/dQ, let us analyze the 
different contributions to Eq. (^. On the right hand side, the first two terms come from 
V • A and the third comes from A ■ V. If A is a constant field, only the third term 
will be nonzero and the 6'-dependent term will be a global factor that will change the 
intensity of dl/dQ for different 6 values, but not its shape. If A is not a constant field. 
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the three terms will survive and as their relative weights are 9 dependent, changes on 9 
also produce changes in the shape of the {dI/dQ)/u! function. 

Moreover, though the contribution of the first two terms is significant, the greatest 
contribution comes from the third term, as shown in Fig. 6 (dashed line), where the first 
two terms of Eq. (H) have been removed. 

We now compare our results with those of two sets of experimental data, in which 
the angular dependence of the photoemission was displayed. In the case of indium,E3 our 
calculated results for the main peak agree reasonably well with the measured ones. The 
measured peak moves from about 0.84 ujp to about 0.92 ujp, as 9 changes from 45° to 
15°, whereas the corresponding calculated peak position moves from 0.77 Up to 0.92 Up 
for the same angle variation. These measurements were not extended to uj > ujp, so we 
cannot perform any comparisons in this region. For Al, on the other hand, the recent 
experimental data of Barman et aS cover the frequency regions both above and below 
Up. Experimentally, the main peak, bellow Up, exhibits no dispersion when 9 changes, 
unlike the calculated results. For uj > ujp there appear in the measured data two minor 
peaks (called features A and B in Ref. |^, which also occur in our calculated results. The 
arrows in Figs. 6 and 7 denote the positions of the high energy peaks, which disperse in 
the same way as the experimental data (feature A of Ref. ^. The vertical line at about 
u/up = 1.07 denotes the position of the minor peak (for 9 = 50° and 9 = 70°) which 
does not disperse (feature B in Ref. |^). 

In conclusion, we have presented an analytical calculation of the photoemission spec- 
trum and its dependence on the photon angle of incidence. The screening of the pho- 
ton field, including nonlocal effects and excitation of longitudinal modes, was calculated 
within a hydrodynamic model with a realistic surface density profile. Our results indicate 
that the main peak, below Up, need not be just due to mulitopole mode excitation. The 
two sets of experimental data with which we compared our results seem to be contradic- 
tory as regards the main peak dispersion. Thus, further experiments on other materials 
are needed to resolve this question. 
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Appendix A 



Here we give the definitions of tlie coefficients Oj and bf. 

ai{z) = iu'^{uj + i'yoiz)) - iu[ATTno{z) + q^ l3o{z)\quad, 

a2{z) = (?{uj + i'yo{z))quad, 

as{z) = iq[Lu/3o{z) - (?{ijj + i'^Q{z))]quad, 



bi{z) = iuq 

62 (-2) = as{z)quad 



l,M^f^-M^) 



quad, 



b^{z) = —AnnQ{z)LU — (u; + i^Q{z)){c q — lu )quad 



bi{z) = uj 
b^iz) = Po{z)uj 



3 no{z} 



En(z) 



quad, 



5 
The coefficients Ai and Bi defined in Eq. (R^ are given by 



AAz) = - 



UJ 



Asi^z) = -iq 



^ 47mo{z) + q^Po{z) 
uj{uj + i'joiz)) 
ujpo{z) 



c^{uj + i-foiz)) 



-1 



Bi{z) 



-''m^J"'^"^' 



z/S 



B2{z) = -iq 



1 _ _ ^ + ^^o{z) 



UJ f3o{z) 



B^{z 
Bjz 



207r 



m^^''^^ Hz) W \ 



{z) = --l^no{z)e^" . 



2)5710 
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FIGURE CAPTIONS 

Fig. 1 Ground state density for aluminum. The solid line is the self consistent calculation 
taken from Ref.ta and the dashed line is the analytical function given in Eq. ( ]25|) for 
6 = 0.66 a.u.. All figures are calculated for r^ = 2.07 a.u.. 

Fig. 2 (a) Real and (b) imaginary parts of the 2;-component of the vector potential vs. 
z calculated for different values of u; < lOp. The normalization constant A^ is the 
transverse component at z\. b = 0.66 a.u., and 6 = 7r/4 were used. 

Fig. 3 Same as Fig. 2 for cj > Up. 

Fig. 4 Variation of the real and imaginary parts of the coefficients Ci, C2 and Di (see Eqs. 
(P^) and (pSD) with u ioi 6 = 0.2 a.u., and 6 = tt/4. Insets: the same coefficients for 
a structureless surface. 

Fig. 5 Variation of the real part of the longitudinal coefficient C2 with u/up for different 
values oi 6. 9 = 7r/4 was used. 

Fig. 6 Photoemission spectrum as a function of uj/ujp for 6 = 0.66 a.u., e = and 6 = tt/A. 
A partial contribution to the spectrum is included as a dashed curve (see text). The 
arrow is at Up/ cos{6). 

Fig. 7 Variation of the photoemission spectrum with 6 for 6 = 1.2 a.u. and e = 0. The 
arrows are at Up/ cos{9). 
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FIG. 1. Ground state density for aluminum. The solid line is the self consistent calculation 
taken from Ref. |l^ and the dashed line is the analytical function given in Eq. (25) for 6 = 0.66 
a.u.. All figures are calculated for r^ = 2.07 a.u.. 



21 



2.00 




r 4 

z [a.u.] 



20 



FIG. 2. (a) Real and (b) imaginary parts of the z-component of the vector potential vs. 
z calculated for different values of lj < ujp. The normalization constant A^ is the transverse 
component at zi. 5 = 0.66 a.u., and 9 = it /A were used. 
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FIG. 3. Same as Fig. 2 for u; > LOp. 
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FIG. 4. Variation of the real and imaginary parts of the coefficients Ci, C2 and Di (see 
Eqs. (ISq ) and (^)) with uj for 5 = 0.2 a.n., and 9 = 7r/4. Insets: the same coefficients for a 
structureless surface. 
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FIG. 5. Variation of the real part of the longitudinal coefficient C2 with uj/ujp for different 
values of (5. 9 = 7r/4 was used. 
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FIG. 6. Photoemission spectrum as a function oiLo/tOp for 5 = 0.66 a.u., e = and = 7r/4. 
A partial contribution to the spectrum is included as a dashed curve (see text). The arrow is 
at LOp/ cos{6). 
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FIG. 7. Variation of the photoemission spectrum with 9 for 6 = 1.2 a.u. and e = 0. The 
arrows are at tOp/ cos{9). 



27 



